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We present the tensor computer algebra package xPert for fast construction and manipulation 
of the equations of metric perturbation theory, around arbitrary backgrounds. It is based on the 
combination of explicit combinatorial formulas for the n-th order perturbation of curvature tensors 
and their gauge changes, and the use of highly efficient techniques of index canonicalization, provided 
by the underlying tensor system xAct, for Mathematica. We give examples of use and show the 
efficiency of the system with timings plots: it is possible to handle orders n = 4 or n = 5 within 
seconds, or reach n = 10 with timings below 1 hour. 
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O I- INTRODUCTION 

'^ , Linear perturbation theory is one of the most successful techniques to extract information from the equations of 

|J^ ' general relativity. It allows us to study the evolution of small deviations with respect to a known background solution 
of the equations, showing the stability properties of that background or how it interacts with weak surrounding fields. 
Linearization is, however, the main limitation of the approach, because we loose the intrinsic nonlinear character of 
general relativity, and hence it is not possible to model the self-interaction or mutual interaction among the linear 

^ i' modes, and other interesting phenomena. Obviously, this could be cured by adding a second perturbative order to 
the analysis, but the large size of the resulting perturbative equations has generally prevented practical computations 

^ until very recently. 

^^ Already in 1980, Cunningham, Price and Moncrief studied and simulated second-order perturbations of collapse of a 

ball of dust (Oppenheimer-Snyder collapse), including the matching to perturbations of an exterior vacuum dj. Since 

^ ' then, other backgrounds have been perturbed at second order. The oscillation modes of a rotating fluid star have been 

^ [ modelled at second order, introducing rotation at linear level [2, l3|. Teukolsky's formalism ^^J for linear perturbations 
'^ ■ of a Kerr black hole was generalized to second order in Q . The gravitational radiation emitted in a close limit head-on 
collision of two black holes has also been studied at second order, describing the final hole as a first-order quadrupolar 
perturbation of the Schwarzschild spacetime [g, [7|. Recently, also perturbing Schwarzschild, second-order quasinormal 
modes have been defined |8|, l9|. Second-order perturbations have been considered in several cosmological scenarios, 
often taking Friedmann- Robertson- Walker as background spacetime; see for example [ifl, [ill, [l2l ^'^'^ references therein. 
We can also find high-order perturbations in various areas of field theory, like the problem of renormalization of a 
r^ ' Lagrangian [ij], though in this context the formalism used is slightly different. We shall later see how the general 
scheme can be adapted to include this as a particular case. 

In references [ij, |l5| we have constructed the arbitrary-order generalization of the Gerlach and Sengupta formalism 
for gauge-invariant and covariant first-order perturbations of any spherical spacetime [16], for any matter model. This 

Vh is the most convenient framework to study nearly-spherical astrophysical scenarios, like collapse of slowly-rotating 

5r , stars, and some types of supernova explosions. The large size of the equations involved forced us to develop specialized 
tools of tensor computer algebra, which were summarized in an appendix of [14| . We have now improved these tools, 
extending their donrain of applicability, increasing their efficiency, sinrplifying their use and fully docunrenting thenr. 
We present version 1.0 of the package xPert for high-order metric perturbation theory. It is based on the combination 
of adapted mathematical algorithms and powerful techniques of tensor computer algebra. On the one hand, we use 
explicit precomputed formulas for the general terms of the perturbative expansions of the most important curvature 
tensors in differential geometry, avoiding slow recursive computations. On the other hand, xPert uses very efficient 
canonicalization algorithms for tensor computer algebra in Mathematica. This allows us to perform computations up 
to perturbative order 10 in several minutes with small computer resources, as we will see in some examples. 

The article is organized as follows. Section [Til reviews the basic concepts of high-order metric perturbation theory 
and presents closed formulas to obtain the n-th perturbation of the most relevant geometric tensors. Section IIIll 
explains how these formulas have been implemented in Mathematica and lists the main commands of the package 
xPert. Finally, we show the high efficiency of the system through some timing examples in section IIVI Section |V] 
contains our conclusions. 
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II. HIGH-ORDER METRIC PERTURBATION THEORY 
A. Basic concepts 

Let us assume a family of A^-dimensional manifolds Al(e) equipped with a metric field 'g{e) depending on a dimen- 
sionless parameter e. This dependence will be taken to be smooth, so that all derivatives with respect to e are well 
defined. Members with parameter s — will be referred to as background members and will be denoted without bar, 
that is A1 = A4{0),g = 5(0). The Levi-Civita connection of g will be called V, and will be denoted with a semicolon 
in indexed expressions. 

The goal in perturbation theory is comparing a tensor field T(e) with its background counterpart T. The former 
is defined on the manifold A1(£) and the latter on the background manifold A4 and therefore we have to introduce 
a family of point-to-point identification mappings (or gauge) 4>e between them to compare. The pull-back (j)*T{e) to 
the background can be now expanded as 

00 „ 



Tie)^T + J2-T^m^ (1) 



n—l 



with all terms in the equation being defined on the background manifold. We have introduced the perturbative 
operator A^ in the gauge 0^, which is actually a derivative: 



Ajm.""*;^'^' 



de" 



(2) 

e=0 



From now on we shall omit the explicit dependence of perturbations on the chosen gauge 0^, to alleviate the notation, 
e-dependent (barred) objects will be also assumed to be pulled-back to the background manifold. The theory and 
formulas for perturbations of changes of gauge have been fully developed for tensors by Bruni and collaborators [I4I , 
and extended in [18] to perturbations of distribution functions in phase space. 

An important consequence of the derivative character of A is that it obeys the Leibnitz rule, whose n-th order 
generalization on the product of m tensors is 

A" [Ti . . . T„] = J2 .,'''. A '' \T,]... A'^'" [r,„] , (3) 

{kij 

where the sum extends to all sorted partitions of n in ?n nonnegative integers (that is, including zero) obeying 
ki + ... + km — n. It is possible to generalize this formula to the case A"[F(S'i, ..., Sm)] for an arbitrary function F of 
m scalar arguments Si, ..., Sm using the Faa di Bruno formula ^.ly| (the n-th order chain rule), generalized to multiple 
arguments. 

We want to do metric perturbation theory, meaning that our fundamental object to be perturbed is the metric 
field gab (abstract indices in any of the tangent spaces of the A4 manifolds will be lowercase Latin indices). It is 
convenient to name the n-th order perturbation of the metric as ^"^^hab = ^"[gabjj with a left-superindex denoting 
the perturbative order. In this notation, the perturbed metric is given by 

5a6(e) ^9ab + J2-( ^"^^^b- (4) 

n=l 

With nonvanishing perturbations of the background metric, it is essential to carefully keep track of index positions. 
For example we have A[ti"] ^ g°-''A[vb] for a tensor field v"-. The notation Av"' could be ambiguous and therefore 
we use square brackets in our notation for perturbations to avoid problems. Should we find the perturbation of a 
tensor with inconvenient index positions, we shall introduce metric factors to correct those positions and then use the 
Leibnitz rule ([3]). 

An important variation of this general perturbative formalism is frequently used in quantum field theory in curved 
backgrounds. It is usually called 'background field method' and decomposes the full metric g as g^b = gab + hab UM- 
Other objects depending nonlinearly on g are later truncated at a given order in powers of h. This can be considered 
as a particular case of the general formalism used here, in which all but the first perturbations of the covariant metric 
are defined to be zero: ^"*/ioh = for n > 2. All formulas in this article can be translated to background field method 
using such a simple restriction. 



B. Perturbations of curvature tensors 

The first step in any perturbative computation in general relativity is expressing the perturbation of one or several 
curvature tensors in terms of perturbations of the metric. This is a straightforward but very long computation for high 
orders, and hence requires the use of computer algebra tools like xPert. There are essentially two ways to perform such 
a computation: either recursively, computing the n-th order perturbation from the previous order, or using explicit 
formulas giving the result directly. xPert implements such explicit formulas for the most important tensors in general 
relativity, as shown in this and the next subsections. Some of these formulas have been already presented in [14| . 

The general term of the expansion of the inverse metric g"'' as a power series in e can be obtained as follows. Let 
us define Hat = Vat ^ 9ab and temporarily drop the indices. Then, the inverse of the metric is given by 



r' 



{9 + H)-\ 



(5) 



which can be expanded around ff = as 
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Reintroducing indices and expanding H itself using ^ we obtain the general term 
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where the sum X](fc ) runs over the 2" ^ sorted partitions of n in m < rt positive integers (not including zero this 
time) such that ki-\- ... + km = n. 

With the general term of the expansion of the metric and its inverse at hand, we can obtain the general term for 
the perturbation of the Christoffel symbols starting from its definition in terms of the metric: 
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where the three-indices tensor perturbation is given by 
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ii^abc — ^ \ 'liable ^ "^ac:b "'bc-.a) : 



(8) 



(9) 



which is symmetric in its last two indices. For n = I formula ([8]) must be read as A[r"bc] — ^^''h°-bc- Note that the 
perturbation of the Christoffel symbols is constructed from the metric perturbations ^"^hab and their first covariant 
derivative at all orders. 

The perturbation of the Riemann tensor of an arbitrary connection (not necessarily a metric connection) is a simple 
application of the Leibnitz rule ([3]) : 
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A"[R,t/]^^b{A"[r\c])-J2 



For a metric-compatible connection we can use 



{k^) 



m 

j2 {'^-m}/i'i<=™...{'=.+i}/i 



fc=i 

getting 



A^-[r«5c]A"-^-[r'^ea]-(a^6). 



If' ••■ i^e-i if'eocb-.a 



(10) 



(11) 



e.+2 -^' *'='/le.e. + ,a <''=- V»^=- . .. '"^'K.^^ '"^'K.^c - (a ^ h). 



s=2 



Other relevant curvature tensors are given in terms of Riemann, and therefore their perturbations can be expressed 
in terms of the latter formula. For instance, perturbations of Ricci are given by 



A"[i?,b] = A"[i?,, 



(12) 



The Ricci scalar R = g°'''Rab can be perturbed using the generahzed Leibnitz rule ([3]): 

A;=0 ^ ^ 



(13) 



and similarly for the Einstein tensor: 

A"[Ga,]^A-[Rat]--J2J2 



n k 



2 A^A^ k\j\(n~j-ky. 



^■>/i,bA'=[/'^]A"-^-'=[i?, 



cd\ 



(14) 



We have not tried to further simplify these two last formulas, as they are already quite efficient and give directly the 
n-th term of the respective perturbative expansions. 

C. Perturbation of the metric determinant 

In Hamiltonian field theory we frequently find the perturbation of the determinant of the metric. This is a basis- 
dependent concept, in the sense that what we are computing is the determinant of the components of the metric in 
a given basis and the result depends on the basis we have chosen. Under a change of basis the determinant changes 
with a squared Jacobian and it is, hence, a density of weight +2. 

The determinant of the metric Qab can be defined as 



det(5Qb 
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(15) 



using the upper antisymmetric density 7y''i-'»« (the overtilde denotes weight +1), whose components in the chosen 
basis are +1, —1 or 0. This object, as well as its lower counterpart ryai...ajv of weight —1, stays invariant under the A 
perturbation. Therefore, the Leibnitz rule ([3]) implies 



A"[det(5,6)] = -l^-^--" ^"-■"^ ^ 



{fc,} 



Finally, this can be simplified using the well known relation 
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We conclude that the n-th order perturbation of the determinant of the metric is always the product of the determinant 
itself times a scalar formed by contraction of metric perturbations. It is interesting to note that such a scalar factor 
contains the product of at most N metric perturbations, and not n, as we might have anticipated by inspection of the 
formulas in the previous section. This is actually the only place in this article in which the dimension of the manifold 
being perturbed plays a role. 

With this formula at hand it is straightforward to give the n-th order perturbation of the volume form ea^ an = 
\det{g)\'/^f]a,...a.. 



D. Perturbations of derivatives 



There are several types of derivatives which may appear in a perturbative computation. We shall study here three 
of those types: partial, covariant and Lie derivatives. 

Partial derivatives are associated to coordinate systems and hence do not change under perturbations. Therefore, 
by construction, they commute with the A operator: A"[Ta] = A"[r]_a for any tensor field T of any rank. 

General covariant derivatives can be perturbed. For instance, the Levi-Civita connection of a metric will change 
when its metric is perturbed. The question arises then about what is the perturbation of the covariant derivative 



of a tensor. Transforming to partial derivatives and ChristofFel symbols, perturbing and coming back to covariant 
derivatives, we get, for an arbitrary tensor density T of weight w, 

A"[r'^^-"'"6i...b,;c] = A"[r'^--"6,...bj., (18) 

I \ ^ f \'n\-\^ai rriai...d...ajn 1 "pii a n rrriai ...d...a„n 1 1 

+ Z^\A [i cdJ f,i...6;,J-i cdA [J bi...b„\i 

1=1 

(A [i cfcjJ f,i...d...6j -i cfcjA [J bi...d...bk\} 

which can be rewritten applying the Leibnitz rule ([3]) as 

A"[r'^^-""6i...b.;c] = A"[T'^--"b,...bJ,, (19) 

+ E ( I ) {E^'[r"'cd]A"-'[r'^-''--™,,...,j-^A'[r%]A"-'[T"--'",,...d...b 



;=i ^ '' i=i j=i 



Finally, the perturbation formula for the Lie derivative along the vector field u of a tensor T of any rank can be 
computed by intermediate transformation to partial or covariant derivatives: 

A"[/:.T] = E ( fc ) C^^^A'^-'m. (20) 

fc=0 ^ ^ 

This expression bears an obvious similarity with the Leibnitz rule, reflecting the fact that both v and T are being 
perturbed but not the Lie structure itself, which is directly given by the differential structure of the manifold and 
thus remains unperturbed. 

III. THE xPert PACKAGE 

The general terms for the perturbation expansions in the previous section are combinatorial in nature and hence 
contain a number of terms which grows exponentially with the perturbative order n. Only through the use of a 
specialized tensor computer algebra system it can be possible to handle them for n beyond 1 or 2. Our work in 
second-order perturbation theory around spherical spacetimes |14l . 11511 motivated us to construct such a system, a 
very preliminary version of which was presented in an appendix of [14| . 

Here we introduce version 1.0 of the free-software package xPert, which implements all those expansion formulas 
and a number of other efhcient tools to work with metric perturbation theory in the standard scenarios of general 
relativity, cosmology, quantum gravity or string theory. This section describes the main commands and features of the 
package, simultaneously constructing a very simple example session. (The In[*J:= and Out[*]:= prompts represent 
respectively input and output lines in Mathematica. Code lines without prompt indicate internal definitions in the 
package.) Next section will provide timings with harder examples. 

xPert is a module of the framework xAct [20,] for efficient tensor computer algebra in Mathematica. xAct is based on 
the use of fast algorithms of computational group theory, which allow very fast canonicalization of indexed expressions 
with arbitrary symmetries |2ll . |22| . Its core tensor packages are xTensor for abstract computations and xCoba for 
component computations, but xAct has also modules for manipulation of the Riemann tensor (Invar [23]) or tensor 
spherical harmonics [Harmonics [20]), and others are under development. These packages share a consistent notation, 
so that they can be run together, and are all free software. 

Following Sectionini which is valid for any background metric, xPert has been developed at the abstract computation 
level, to avoid dealing with choices of coordinate systems or frames. Therefore, it runs on xTensor, using its notations 
and conventions. A tensor field TJ' is denoted by T[-a, b] and the partial derivative operator da is represented as 
PD[-a] , such that daV^ is PD[-a] [v[b]] in the system. Let us start by loading xTensor. 

In[l] := «xAct' xTensor' 

(Version and copyright messages) 



We first define our background structure: a four-dimensional manifold M whose tangent vector space will have 
(abstract) indices {a,b,c,d,e}, 

In[2] :^ DefManifold[ M, 4, {a,b,c,d,e} ] 

Then we define a metric tensor field g with negative determinant and associated Levi-Civita covariant derivative CD, 

In[3] :^ DefMetric[ -1, g[-a,-b], CD, {";","V"} ] 

(Info messages on construction of associated tensors) 

We have provided the symbols {";"," V " } to format the derivative in postfix or prefix output notation, respec- 
tively. DefMetric automatically defines all tensors normally associated to the metric or its connection, such as 
Christof f elCD [a,-b,-c] , RiemannCD [-a,-b,-c,-d] , EinsteinCD [-a,-b] , Detg[] , and so on, with obvious mean- 
ings. We can define other tensors with the syntax 

In[4j := DefTensor[ MaxwellF [a,b] , M, Antisymmetric [{a, b}] , PrintAs->"F" ] 

The arrow -> is the Mathematica representation for an optional named argument. 

Now we load xPert (this would also load automatically xTensor if it wasn't already in memory): 

In[5j := «xAct ' xPert ' 



Package xAct'xPert' version 1.0.0, {2008, 6, 30} 
Copyright (C) 2005-2008 David Brizuela, Jose M. Martin-Garcia 
and Guillermo A. Mena Marugan, under GPL 



This adds several new commands and reserved words to the system, of which we shall here describe the four most 
important, namely Def MetricPerturbation, Perturbation, ExpandPerturbation and GaugeChange. 

A perturbative structure having metric g as background and the tensor h as its perturbation is defined using 

In[6] := Def MetricPerturbation [ g, h, £ ] 

which also identifies e as the perturbative parameter of the expansions. From now on, the rt-th perturbation of the 
metric g[-a,-b] will be denoted as h[LI [n] ,-a,-b] , where LI is the iTensor head to denote so-called 'label indices', 
that is, indices with no vector space association. Labels can be considered as general non-geometric purpose indices. 
The A operator of section |II] is represented by the head Perturbation. It has two arguments: the background 
expression being perturbed and the perturbative order (with default value 1): 

In[7] :^ Perturbation[ MaxwellF [a, b] , 3 ] 

Out[7j ;= A^lF"''] 

Note that the tensor is represented with its symbol F and that the perturbation order is an exponent of A, as in 
section |TT1 Following normal Mathematica, the output is formatted for the sake of clarity, but the internal notation 
is still the same. Perturbation acts mainly as a wrapper for tensor expressions, but has been instructed to evaluate 
them under certain conditions. First, it automatically combines perturbative orders of composed heads (symbols with 
an underscore are named patterns in Mathematica): 

Perturbation [ expr_, ] := expr 

Perturbation[ Perturbation[ expr_, n_ ] , m_ ] := Perturbation[ expr, n+m ] 

Being a derivative, Perturbation is linear and gives zero on the delta identity tensor and constants: 

Perturbation [ x_ + y_, n_ ] := Perturbation [ x, n ] + Perturbation [ y, n ] 



Perturbation [ delta [a_, b_] , n_ ] := 

Perturbation [ expr_?ConstantQ, n_ ] := 

The Leibnitz rule is also automatic, and has been implemented following equation ([3]) for any number of factors and 
any perturbative order, using fast algorithms to compute partitions implemented in xPert. Perturbation commutes 
with partial derivatives of general expressions and with any covariant derivative of a scalar expression: 

Perturbation[ PD[-a_] [ expr_ ] , n_ ] := PD[-a] [ Perturbation[ expr, n ] ] 

Perturbation [ CD_?CovDQ [-aj [ expr_?ScalarQ ] , n_ ] : = 

CD [-a] [ Perturbation [ expr, n ] ] 

The index of the derivatives is required to be always covariant, to avoid a metric mismatch, and that is implemented 
through a pattern index -a_ . Finally, Perturbation does not change the density weight of the perturbed expression: 

WeightOf [ Perturbation [ expr_, n_ ] ] := WeightOf [ expr ] 

The Def MetricPerturbation in In[6] defines special rules for the metric g and its perturbations h with covariant 
indices: 

Perturbation [ g[-a_,-b_], n_ ] := h[LI [n] ,-a,-b] 

Perturbation[ h[LI [n_] ,-a_,-b_] , m_ ] := h[LI [n+m] ,-a,-b] 

With the setup and internal definitions so far we can now perform computations like this second order perturbation 

In[8] := Perturbation[ g[-a,-b] RicciCD[c,d] + RiemannCD [-a,-b, c,d] , 2 ] 

Out [8] ;= 2h\bA[R,d]+9abA^[Rcd] + A^[Rab^''] + h\bRcd 

Actually, we could now proceed to perform any computation in metric perturbation theory by decomposing the 
curvature tensors in partial derivatives of the metric and using the code already given recursively. Only the definition 
A [5"''] = — {i}/ja'> would be missing. However, that would be highly inefhcient already for moderate perturbative 
order n. Instead, we shall use the expansion formulas of section HIl which allow the nonrecursive construction of 
perturbations at any order n. 

Formulas P^ [^D]) for derivative expansions and formulas Q [51 [Til [12 [HI [HI [IS|) for the relevant curvature 
tensors have all been encoded in a single command called ExpandPerturbation, the most powerful part of xPert. 
ExpsindPerturbation takes any expression and replaces the arbitrary-order perturbations of known background ob- 
jects by their expansions in terms of metric perturbations, but only if those objects have their indices in the appropriate 
positions. For example the perturbation of the Einstein tensor has only been stored for covariant indices. In all other 
cases there is an internal call to the xTensor function SeparateMetric, which introduces metric factors to bring all 
indices to their natural positions, which arc those given at definition time. To show how this works, we perform an 
explicit metric separation by hand (symbol 7o represents the previous output): 

In[9] :— Perturbation[ EinsteinCD [a,b] ] 

Out[9] := A[G"'''] 

In[10] := SeparateMetric [ ] [ 7, ] 

Out[10] := Gcd g""" A[g"^] + g''^ ( g'"^ A[Gcd] + G^d A[g^''] ) 

Now ExpandPerturbation can expand the perturbation of the Einstein tensor with covariant indices, and the per- 
turbation of the inverse metric. 

In[ll] := ContractMetric [ ExpandPerturbation [ 7„ ] ] 

Out[nj ;= -Gc''h^'"'- - G^'ch^^" + \g''^h^'"^Rcd ~ ^h^''''R - i/i^^c''''" 

1 ulcb ;a I 1 J-,lb -^c-.a i 1 Llc6;a i 1 ulca.b 1 /jl6a;c 

2 "' ;c ~r 2''' c ~r 2''' ;c ~r 2"' ;c 2 1*^ 



AbsoluteTimlng [ 
Perturbation [EinsteinCD [-a, -b] , 2] // ExpandPerturbation // ContractMetric // 
ToCanonlcal ] 



{0.783872 Second, h\^ h^ -^ R,^ + ^ g,^ h^ -^ r^^ _ g^^ h^/ h^ '^'^ R^, - ^ h\^ R + h^ -^ h^ 



cd; b; a 
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„ r.^lc;di,^le i-^lcdi,^le vlcdi,^l ;e ,,,, h^ h^ '^'^' ^ n h^ Vi^ '^'^'^ 1 

^ab "■ c '^ d ;e + ^ab "^ "^ c ;d;e ~ ^ab '^ ^1 cd ; e + "g" ^Jab "■ ce;d ^^ " ^ ^ab " cd; e '^ J 

FIG. 1: The second-order perturbation of the Einstein tensor is constructed and canonicalized in less than one second. The 
blue labels of the h tensors denote the perturbative order. 

I l^,abuld :c I 1 ^ab Icicle 1 ^abl,lc \d l^abuldc 1 l^abulc :d 

+ 49 " c ;c + 4ff n, ;d;c - 45 n c' ;c — 2^ " ;c;d + 45 "■ c ;rf 

The xTensor command ContractMetric has been used to absorb all possible metric factors. Finally, ToCanonical 
moves indices around bringing equal terms together, 

In[12] := ToCanonical [ 7, ] 

Out[12] ;= -G^'^h^'^c - C^h^K + y^^h^'^'^Rcd - ^h^'^^R - ^h^^^"'''' 

+ i/i"='";c + i/i^°'='^c - ^h^'''''''';c - y''''h^'"^-c-d + \g''^h^''c''^;d 

Figure [T] shows the output of the second perturbation of the covariant Einstein tensor, computed with the same 
combination of commands. 

Another useful command implemented in xPert is GaugeChange. A e-dependent family of diffeomorphisms on a 
manifold can be parametrized by an infinite collection of vector fields f"'^°, one per perturbative order, and the action 
of these on any background tensor at any order has been given by Bruni and collaborators [131 . The general case has 
been implemented here and can be used as follows. First define the family of generator vector fields on the manifold 
M transforming from the current gauge to a new gauge: 

In[13] :^ Def Tensor [ ^[LI[n],a], M ] 

The third-order perturbation of F°-^ can be changed to the new gauge using 

In[14-J :— GaugeChange [ Perturbation[ MaxwellF [a,b] , 3 ], ^ ] 

Out[14j := A^[F'''']+3C^iA^[F'''']+3C^iC^iA[F'''']+C^iC^iC^iF'''' 
+3 C^i C^2 F'''' + 3 £^2 A [F''''] + C^3 F"^ 

We finish this section by coming back to the problem of perturbation theory using the background field method, in 
which all but the first metric perturbations vanish, as stated in section [TTl This can be easily implemented setting 

In[15] :^ h[LI[nJ ,a_,b_] :=0 /; n > 1 
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FIG. 2: Canonicalization timings (in seconds) for the perturbation of the Riemann tensor at perturbative orders n = 1...10 
(lower, red line). Also shown number of terms in the expression (upper, blue dashed line). Both lines are clear exponentials in 
n, with the timings growing slightly faster because terms with larger n are harder to canonicalize due to their larger average 
number of indices. 

IV. EXAMPLES 

This section presents several examples of computations, focusing on the dependence of their timings on the per- 
turbative order and the number of objects being perturbed. The intrinsic combinatorial nature of the problem will 
always imply exponential dependence, but we will see that the timings in xPert are short enough to handle all useful 
cases. These examples have been performed using a Linux box with a 3 GHz Pentium IV processor and 2 Gb of RAM. 

In perturbation theory the overall level of complexity is mainly determined by the perturbative order n. It affects 
the computation in two different ways: on the one hand the expressions to manipulate are sums with a number of terms 
which grows exponentially with n; on the other hand each term is a product of objects and the number of factors also 
grows (typically linearly) with n. Canonicalizing a sum of terms is obviously a linear process because each term can 
be dealt with independently (this is an ideal scenario for parallelization), but canonicalization of a product of objects 
is naturally factorial in the number of indices and this could prevent any practical computation. The canonicalization 
algorithms in xTensor are fast enough to render the problem effectively polynomial in the number of indices, allowing 
us to deal with expressions of a few dozen indices in hundredths of a second. Figure [5] shows the number of terms 
and the timing of canonicalization of the perturbation of the Riemann tensor at different perturbative orders. (The 
timing of construction of the expression is negligible in comparison.) The 10-th order perturbation is canonicalized 
in slightly less than 20 minutes and contains 44544 terms. We see clearly the exponential growth of both curves, but 
with manageable timings. Third-order perturbation expressions can be manipulated in 1 second. 

Other possible sources of complexity in perturbative computations are the expansion of perturbations of a product 
of tensors and the expansion of perturbations of a function of a number of scalar arguments. Our implementations 
of the nth-order Leibnitz rule and Faa di Bruno formula, respectively, are fast enough to neglect their timings in 
comparison with those of canonicalization. Figures [3] and [3] display example timings for those problems. The Leibnitz 
rule is simpler than the Faa di Bruno formula and produces faster results, also taking less memory. 

Overall, we see that we are limited in size by RAM memory, which allows us to work with up to 10^ terms roughly 
with a few Gbytes, corresponding to n = 10 approximately. Time limitations come mainly from the canonicalization 
process (other expansions are faster): Within seconds we can manipulate all equations up to orders n = 4 or n = 5. 
The n = 10 equations require canonicalization times of the order of 1 hour. This gives an idea of the power and 
efficiency of xPert, and what can be achieved with it. 



CONCLUSIONS 



Linear perturbation theory is an important tool in general relativity, nowadays playing a complementary role to 
numerical relativity [2J| . Its precision and applicability can be extended by adding higher-order perturbations in the 
computation, transforming the intrinsic non-linearity of the theory into sources for the linearized equations of motion. 
In most cases these sources are huge, and hence practical calculations require the use of specialized computer algebra 
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FIG. 3: Timings of expansion (in seconds; red lines, stars) and number of terms (blue dashed lines, diamonds) of the perturbation 
of the product of m factors, for different perturbative orders n = 1...10 [formula ©J. Different lines correspond to increasing 
values of m, from m = 2 to m = 10 starting from below. All practical cases stay below one second. 
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FIG. 4: Timings of expansion (in seconds; red lines, stars) and number of terms (blue dashed lines, diamonds) of the perturbation 
of a scalar function of m scalar arguments, for different perturbative orders n — 1...10. Different lines correspond to increasing 
values of m, from m = 1 to m = 10 starting from below. We include only those cases which can be handled with 2Gb of RAM 
memory, corresponding to a few tens thousand terms. For example, for n = 10 we can only handle up to m = 4. 



to construct and manipulate them. 

In this article we have introduced version 1.0 of the package xPert for high-order metric perturbation theory. 
This is a very hard problem, of intrinsic combinatorial nature, and hence exponential in the perturbative order n. 
We have developed and coded up highly efficient algorithms and tools, based on the combination of explicit non- 
recursive formulas to build the equations at any order n, and powerful algorithms of computational group theory to 
manipulate them. xPert can work both with a perturbative scheme in which the metric is expanded as an infinite 
series of perturbations (the usual approach in general relativity) and with the alternative scheme in which the metric 
is written as a background plus a single 'deviation' term (as used in quantum field theory). xPert allows working with 
orders n = 4 and n = 5 in seconds, and we have shown that it is possible to reach n = 10 with average timings below 
1 hour with moderate computer resources. The main limitation (both in time and memory) is the huge number of 
terms produced, but this can be easily dealt with using parallelization, because different terms in the expression can 
be treated independently by different processors. 

xPert is a module of the general framework xAct [201] for tensor computer algebra in Mathematica, sharing a 
consistent notation with the whole system. This means that once the perturbative equations have been constructed, 
there are several hundred additional commands available for tensor manipulations, including particularization to a 
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special background metric, and several thousand more for generic manipulation in Mathematica. In particular, if one 
wants to convert the produced equations into C or FORTRAN code, it can be easily achieved with the CForm and 
FortrsmForm commands in Mathematica. 

xPert has been already used in several research projects: Our construction of a general formalism for high-order 
metric perturbations around spherical backgrounds [14| and its gauge-invariant form [15]. In a cosmological setting, 
high-order perturbations of non-linear radiation transfer have been studied using xPert, both in kinetic theory [l8| 
and in a fluid approximation [25| , as well as first-order perturbations of scalar field inflation for anisotropic spacetimes 

M- 

The package xPert is free software distributed under the GNU general public license, and can be download from 
the webpage http: //metric, iem. csic. es/Martin-Garcia/xAct/xPert/| . 
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